QCD SUM RULES AND THE DETERMINATION OF 
LEADING TWIST NON-SINGLET OPERATOR 
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Abstract 

We use QCD sum rules to determine directly the leading-twist non-singlet operator 
matrix elements based on calculations of three-point correlator functions in configuration 
space. We find a different result from that obtained by integrating the structure functions' 
expressions obtained by Belyaev and loffe based on calculations of four-point correlators 
in momentum space. The origin of this discrepancy remains unclear. 

1 Introduction 

The determination of the functional dependence of the nuclear structure functions on the 
Bjorken scaling variable x via QCD sum rules has been pioneered by Belyaev and loffe |jl|, 0]. 
They considered the four-point correlator 



d'xdSd'z e"^^e'^^y-^\T{r^{y)j,{x)j,{Q)n{z)})o 



(1) 



where rj is the standard three-quark current with proton quantum numbers and is the 
electromagnetic current (see Fig. |I]). The sum rules they obtained have a limited applicability 
region restricted to intermediate x far from the kinematical boundaries x = and 1, and should 
be used for not very large ^ 10 GeV^. 

In this paper we propose another method to calculate the quark distributions by evaluating 
their moments. More precisely, we shall use QCD sum rules to determine the leading twist-two 
non-singlet operator matrix elements (OMEs) contributing to deep inelastic scattering from nu- 
cleon targets which, by the moments sum rules, are related to the quark distribution moments. 




Figure 1: The lowest-order contribution to the four-point correlator of Equation |l]. 




Figure 2: Applying the OPE for the upper propagator in the four-point correlator approach 
leads to the three-point correlator approach. 

In principle, this method should be equivalent to that of since it is based on calculation of 
three-point correlation functions which can be produced by performing the OPE on the two 
electromagnetic currents in the four-point correlator approach (see Fig. However, it does 
not have the problem as it is calculated perturbatively from the perturbative dependence 
of the coefficient functions and the anomalous dimension of the operators in the OPE series. 
It is crucial to note that, assuming the validity of the conventional OPE as was the case in 
the calculations for higher moments do not reveal logarithmic divergences. Thus it is not 
necessary to use the external field method contrary to the case of the second moments where 
the logarithmic divergences are present 0. [] 



2 The method 



Following the standard procedure of QCD sum rules, we calculate the three-point correlation 
function 

T,,(p) = |d^(x)e*^^'(0|T{r7,(x)6^„,...«„r7,(0)}|0) (2) 

where ?7i(x) is the interpolating current for a proton originally suggested in [^]. 

= (nr(x)C7;.n,(x))(757^4(x)),e'^'''= (3) 

C is the charge conjugation matrix, superscript T means transpose, i and j are spinor indices 
while a, h and c are colour indices. As to the operator 0/3a^...Q„, we choose the twist -two 
non-singlet composite operator given by 

= V-S ^dS'^{y)^pD^,...D^J>{y) (4) 

where Dq, is the covariant derivative 

^.\y = d.-ig{^)Al{y) = a„-^|(y)?/,G;jO) (5) 

^In fact, performing the calculations in the case of the second moment would lead to integrals of the form 
/ d'^xe^P''^-^^ where m > 2. These integrals when Borel transformed give a divergent contribution. 
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in Schwinger gauge where yaA^{y) = 0, and S denotes the symmetrization over the indices 
(3, ai, . . . , an with subtraction of traces in order to make the operator O of definite spin. 

The operators Of3ai...a„ form an irreducible representation of the Lorentz group. Their 
reduced matrix elements are defined by 

{N{p)\dpai...ajN{p)) = An p^p"'. . .p"'" + terms containing g^^^^ , (6) 

and our aim is to estimate these OMEs' A"- by the QCD sum rules method. 

First, we consider the "phenomenological" representation of the three-point correlator which 
is expressed in terms of physical intermediate states as 

{OlrjilN^p)) X nucleon propagator X { N{p)\0f3ai...an\^{p)) 

X nucleon propagator X { N{p)\f]j\0) (7) 

where nucleon propaqator = i f^™, and m is the mass of the nucleon A^. In addition to this 
double pole term, there are also single pole (continuum) terms suppressed relative to the former 
with power of the Borel parameter (exponentially). One can check that the combination 
of invariant functions which enters the coefficient at the structure pp^p""^ . . .p"" or p^p'^^. . .p"" 
in Equation |^ coincides (up to a numerical factor) with the combination of invariant functions 
at the structure in the spin-averaged matrix element { N{p)\0i3ai...an\N{p)) i-e. 

with the reduced matrix element we are interested in, and we get: 

Phen = —T. :rT^2An(p + rn)p^p°'^ . . .p"" + other structures, (8) 

(p"' — m^)'' 

where A is the constant defined by 

mNip)) = Xv, (9) 

and Vp is the proton spinor satisfying {p — m)v.p = and the normalisation J2poiarisations r ^^'^^ = 
2m. 

The different tensor structures emerging from the double pole term can be used to construct 
a sum rule and extract An, but in our case the sum rule from the structure pp^p"'^ . . .p"'" 
is preferred. Firstly, this is because it contains the maximum number of momenta in the 
numerator and thus improves the convergence of the OPE series and diminishes the background 
contribution of excited hadronic states compared to the lowest state (proton) contribution of 
interest Secondly, the structure pp^p°'^. . .p"" conserves chirality: this is also a merit since 
for structures conserving chirality one can calculate a larger number of terms in the OPE series. 
Finally, for this tensor structure, several simplification tricks proved to be possible which made 
the calculations for the "theoretical" part manageable. 

We consider now the "theoretical" part of the sum rules. In the process of the calculations 
we shall take into account only the operators of dimension d <6 and calculate only the invariant 
functions at the structure ppf^p°'^ . . .p"" . The diagrams corresponding to this tensor structure 
are depicted in Fig. ^ where the unit operator contribution corresponds to Fig. the gluon 
condensate {G^^^C^^^) to Fig. ^-1 and the quark condensate {i^i')^ to Fig. ^|m. For massless 
quarks, it is convenient to work in configuration space and in order to calculate the gluonic 
condensate contribution it is far easier to work in the fixed-point gauge for the gluonic field 
z^A;{z) = 0. 

The standard choice of the reference frame is to identify the origin {z = 0) with the coor- 
dinate of one of the currents, so we have the following two choices. 
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Figure 3: The "theoretical" sides of the sum rule for the three-point correlation function. The 
wavy lines denote the current, and the crossed circles denote the operator insertion. The solid 
and dashed lines denote quarks and gluons respectively. 



1. The operator O is inserted at the origin (Fig. Ha): The advantage of this gauge is that, 
because covariant derivatives may be replaced by ordinary derivatives at the origin, the 
graphs where gluons originate from the operator (Fig. 0i, j, k, 1) do not contribute. This 
simplification is, however, illusory and one should be very careful in doing the calculations 
in this gauge choice. In fact, because the operator contains derivatives, it is impossible to 
put {z = 0) from the very beginning since the quark loop is determined now by derivatives 
of the type {^)S{z,y) where S{z,y) denotes the propagator between the two points y 



and z and is given by [§, |^, §] : 



Sfj'^{z,y) 



27r2 {z - 
27r2 2 



yY ^ 167r2 2^'''^\2 
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[z y 



yf 



~y) 



{z-yY 



(10) 



where A,B are colour indices, i,j are spinor indices, 
Gell-Mann matrices. 



and A" are the 



Since one has to differentiate first with respect to z and only put 2; = afterwards, we 
can see that the "non-translation-invariant" terms in S{z, y) would give non-zero results 
leading to new integrals to be done. 

2. The current fj is inserted at the origin (Fig. ^): Here the derivatives are applied at y and 
no derivatives in the current 77 at so we can put z = Q from the beginning. Looking 
again at the expression for S{z, y) we see that the graphs in Fig. ^f and h contribute 
nothing but those of Fig. j, k and 1 do contribute since the covariant derivative D 
contains operators of the gluon field (c.f. Equation |^). 

At first glance, this gives rise to some doubts as one might suspect that the gluons 
originating from the operator (e.g. Fig. ^) correspond to gluons emitted from the upper 
propagator in the four-point correlation function formalism (Fig. §a). If this were true 
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(a) (b) 

Figure 4: Two equivalent choices for the origin. 





(a) 



Figure 5: The contribution of graph (a) is of higher twist while graph (b) may contribute to 
leading twist effects. 



then it would lead to inconsistencies, because in Schwinger gauge the gluon field is 
expressed in terms of the field strength tensor G^u and the standard analysis of OPE in 
DIS shows that this would be of higher twist effect while our operator is of leading twist! 
However, as the choice of gauge in |I] shows, such higher twist operators do not contribute. 
Thus gauge invariance requires that these operator-originated gluons have their origin in 
the right-hand upper propagator (Fig. ^d) and no inconsistencies persist. 

We opt for the second choice of gauge: the origin at the right-hand vertex of the diagrams 
of Fig. 1^, y as the point where the operator 0/3ai...a„ is inserted while x is the point where the 
nucleon is created. For this gauge choice, the calculations are easier and can be done in two 
steps. We do the y-integration first to obtain expressions of "full" propagators in the presence 
of the operator Oj3ai...a„, leaving the a;-integration till the end. 

As an example, doing the ^-integration in Fig. ^ results in a full propagator 
A _B 



(0) 




— L OJ a {/ (x-y)i IP^a-i ■ ■ -Ua^ ^4 

where we have used the relation 

As we are interested only in the structure , the derivatives can be done giving 
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Essentially the same techniques are used to evaluate other "full" propagators in the presence 
of operators even though the calculations are more complicated. After this, one can evaluate 
the graphs: for example, we get for the bare graph where the operator is inserted on a u-quark 
line the result: 

\ 9/ -48-36n . 



In order to get the expression in the momentum representation, one should now perform 
the x-integration using the formula 

^4 ■( _'\\nn'i:~2n 2 




(a;2 - ?e)" T{n - l)T{n) 

For the other graphs, the calculations were performed using the program Mathematica for 
symbolic calculus (we used the FeynCalc package). The results of all the graphs are listed 
below (the letters a,b,. . . correspond to those of Fig. |^, the superscripts u and d denote that 

the operator is acting on u- and d-quarks respectively while the subscripts u and d mean that 
a u- or d-quark line where the operator is not inserted emits a gluon). 

1/1 3 -15 4 \ . . 

C*^'" = ^(l92(^ + WT2))i? 

(tkcr, = ( —1— + —il—] 1. (15) 

ur = ^'<^'>/ 1^1 I I U (16) 

TT^ \A32{n-l) U4n 144(n+l) 432(n + 2) / p2 ^ ^ 

7r4 1^192(72 + 1) 48(71 + 2)^^2 ^'^ 

(or = / I I 5 \ 1 

vr^ I^144n 192(n + 1) 288{n + 2)Jp'^ ^ ' 

irtjT = + ^ + U (19) 

^ 7r4 \^432(n-l)192(n + l) 864(^ + 2)^^2 ^ 

(fc)" = (20) 

^ ^'^ 7r4 1^192(n + l) ^ 192(n + 2)j p2 l^^^ 

(/)" = -5 + M - (22) 

l^576(n + l) 72(n + 2)yj92 ^ > 

(m)" = (23) 
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1 -3 1 



H — -, — - — r H 1 In —p^ 



{hkc 



Att^ \2{n+l) 2(n + 3) n + 

9'{G') ( -1 1 \ 1 

vr^ 1^64(n + l) 6A{n + 2)Jp' 



,2 



^ ^ TT^ \432(r2-l) 192n 288(n + l) 1728(n + 2) 



vr^ \^192(n + 1) 96(n + 2) J 



\^192ri 192(n + 1) 576(n + 2)/ p2 

^ l^432(n-l) 192(n + l) 216(n + 2);j92 

. _ 9'{G') ( -1 1 \ 1 

192(n + l) 192(n + 2)yp2 



{k) 



ilY = ^ ^ ' ^ + 

^' \^576(n + l) 288(n + 2) y p2 

3 p* 

Summing all the graphs we arrive at: 

Thed" = ( + — ^ + + — ) p^ln - p' 

327r4Vn + l n + 2 n + 3 n + Aj^ ^ 

g^(G^) ( -4 -2 



,2 



327r4 \27(n-l) 3(n+l) 27(n + 2) / p2 

Theo" = + ^ + p'lr.-p' 

327r4 Vn + 1 n + 3 n + AJ ^ ^ 

g^{G^) / -4 -2 28 \ 1_ 

^ 327r4 1^27(n- 1) ^ 3(n+ 1) ^ 27(n + 2) J ^ 

-4/ s2 1 

3 p^ 

Equating with Equation |^ and doing the Borel transform, we arrive at 



m? " 2\lj~' i'\n + l ' n + 2 n + 3'n + 4 



+ 



6/41 21 41 



M^V 27n-l 3n+l 27 n + 2 



M2 . ^2/A^2 r, / 1 -3 2 



" m2 2\%'' rVn + 1 ' n + 3 ' n + 4 

6/41 21 28 1 
+ TT7 7-:t^— + 



M^V 27n-l 3n+l 27 n + 2 
\ 



a" 



^ 3M6 j 
where m is the nucleon mass, 

a = -{2nf{Q\ijm , 



6=(27r)2(0KG%,G%,|0) 

TT 

and 

= 2(27r)^A2 . 



3 Comparison with other results 

As mentioned in the introduction, Belyaev and loffe [|1| determined the valence up and down 
quark distributions in the proton by considering the four-point correlator T^^, (Equation |l|) 
where x = ^ is fixed and not too close to the boundaries x = 0, 1 and = —q^ is assumed to 

be only moderately large (Q^ ^ lOGey^) since the logarithmic corrections [asfln -rr — )]"■ were 

^^QCD 

not summed. Assuming the legitimacy of using the conventional OPE for the imaginary part 
of the correlator, they calculated Im T^^ taking perturbative and non-perturbative corrections 
into account, performed the Borel transformation over the parameter keeping only leading 
powers in an expansion in (^) and hence projecting only the leading twist contribution. They 
obtained the following sum rules (the effects of the continuum, the anomalous dimensions and 
the terms oc ag < qq >^ are not shown): 

xu^{x,Q^)+M^A''{x,Q^) = ^e™'/^^'{4a;(l-a;)2(l + 8x) 

xd,ix,Q^) + M^A\x,Q^) = -^e™'/^^' {4x(l - 

2A^ 

H 7 \ X x^ '> (31 

V 27x 6 12 54 ' ' ^ 

The standard analysis of DIS would, in principle, allow one to calculate the twist-two OMEs 
simply by taking moments of the structure functions 

Al = r dxx''u^{x,Q'^) (39) 




Ad 



dxx''d,{x,Q^) (40) 

Jo 



and formally one obtains 



m^"-" 2Xij' l"Vn + l ' n + 2 n + 3 ' n + 4. 

1 71 19 1 97 1 



27n-l 6n 12 n + 1 108n + 2. 




(41) 



m2 " 2Ai lVn + 1 n + 3 n + A 



4 1 71 11 1 -7 1 

+ 77- - 



27n-l 6n 12ra + l 54 n + 2 



(42) 



The term proportional to in Equation ^ corresponds to the 6{l — x) piece which was omitted 
in the structure function expression in Equation |3l apphcable at intermediate values of x. 
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4 Discussion 



In comparing Equations ^ and |36| with Equations ^ and we see that the unit operator and 
< qq >^ contributions agree while there is a discrepancy in the < > term. In Belyaev's opin- 
ion PI, it is not necessary that (|35|, |36|) and (O, B2[) give the same answer because the sum rules 



( p7| , pq ) are used to find the quark distribution for intermediate x only, and one cannot integrate 
them analytically over the whole region x G [0, 1] because there can be singular contributions 
of the OPE near the points x = 0, 1. However, in a formal (rather than phenomeno logical) 
sense, one can compare the two methods. Since we can identify the singularities at x = 0, 1 in 
the approach of we can formally evaluate the moments and compare them with the OPE 
analysis. We have not been able to establish the origin of the discrepancy although, as was 
noted in |jlO|, the phenomenological analysis of the OMEs [jll], |12| changes only slightly using 



the new < > term expression. However the difference is still of importance and we think it 
essential to understand the origin of the discrepancy in the two methods. We hope this paper 
will stimulate further investigation of this puzzle. 

Acknowledgement: I am very grateful to G.G. Ross for help in this analysis. 
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